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CHAPTER  I 

INTRODUCTION  AND  SlT^l'-IARY 

1.   Introduction 

The  steady  state  spatial  distribution  of  neutron  fluxes  in  matter  is  a 
problem  of  prime  Importance  to  the  designer  of  chain  reacting  piles,  and 
has  been  the  subject  of  very  much  theoretical  and  experimental  work.   The 
case  in  which  the  neutron  energy  is  large  compared  to  that  of  the  nuclei 
constituting  the  medium  (the  "slov/ing  down"  case)  is  of  first  importance 
and  is  rather  extensively  investigated.   The  alternate  case  in  which  the 
neutrons  are  in  energy  equilibrium  with  the  medium  is  of  lesser  importance 
and  comparatively  little  experimental  v/ork  has  been  done,  although  the 
fundamental  theory  is  well  in  hand.   The  cross-sections  for  neutron  scat- 
tering and  absorption  appropriate  to  the  first  case  are  the  "ordinary" 
energy-dependent  cross-sections  (or  simple  functions  of  these)  which  have 
been  so  extensively  investigated  by  various  neutron  spectrometry  groups. 
The  cross-sections  appropriate  to  the  second  case  are  more  vaguely  defined. 
Perhaps  as  adeqxiate  a  definition  of  these  parameters  as  any  is  the  "ad  hoc" 
one  -  those  numbers  which,  when  inserted  in  the  mathematical  theory  of 
particle  transport  yield  proper  solutions  for  the  neutron  flux.   It  can 
be  demonstrated  that  these  numbers  are  not  in  general  the  average  values 
of  the  energy  dependent  cross  sections  weighted  by  the  appropriate  energy 
distribution  function  for  the  ambient  temperature. 

The  "ad  hoc"  definition  suggests  a  means  of  experimentally  measuring 
the  numbers.   One  requires  a  medium  -  of  simple  geometry  so  that  the  trans- 
port problem  is  soluble  -  in  which  a  flux  of  thermal  equilibrium  neutrons 


is  maintained  and  in  which  the  flux  can  be  measured  at  various  points  with 
neutron  detectors «   Since  two  parameters  are  to  be  determined,  two  independent 
solutions  for  the  flux  must  be  available  -  but  this  is  easily  provided  for 
by  considering  (l)  a  region  "remote"  from  discontinuities  in  which  the 
solution  is  asymptotic  and  (2)  a  region  adjacent  to  a  discontinuity  -  say 
a  "black"  boundaryo   One  measures  the  flux  at  a  sufficient  number  of  points 
in  each  region  to  define  the  spatial  functions.   The  transport  problem  is 
solved  many  times  for  various  values  of  the  parameters,  and  if  fortune  is 
good  one  pair  of  values  will  be  found  to  yield  solutions  which  fit  the  ex- 
perimental functions  in  both  regions. 

To  particularize  ■=  the  chosen  medium  is  normal  water,  the  geometry 
one°dimensional  consisting  of  a  tank  many  mean  free  paths  in  its  dimensions, 
traversed  by  a  uniform  thermal  neutron  current  arriving  at  one  face  and 
absorbed  by  a  layer  of  cadmiuir;  at  the  opposite  face.  The  neutron  flux  is 
measured  at  points  along  a  line  normal  to  the  cadmium  facing  and  extending 
from  the  facing  well  into  region  (l)o 

The  one-dimensional  transport  problem  is  integrable  in  closed  form  for 
region  (l).   This  solution  is  matched  to  the  experimental  data,  yielding  a 
characteristic  equation  in  say  the  total  cross-section  and  the  scattering 
cross-section  ■=  thus  the  scattering  section  may  be  eliminated.   The  trans- 
port problem  is  then  integrated  numerically  in  region  (2)  for  various 
values  of  total  cross  section  and  that  value  selected  which  yields  best  fit 
with  experimental  data.   Unfortunately,  this  straight-forward  procedure 
could  not  be  carried  out„   It  was  found  that  numerical  integration  produced 
a  function  which  fit  the  experimental  data  astonishingly  well,  but  that  the 
values  of  the  parameter  which  yielded  this  good  fit  were  variable  with 


respect  to  the  integrating  procedure.   Further,  it  was  found  that  each 
improvement  in  the  integrating  plan  led  monotonously  to  a  reduction  in 
the  total  cross  section,  indicating  that  the  process  of  improving  the 
operation  was  not  near  convergence.   The  rapid  increase  in  computing  labor 
made  it  necessary  to  abandon  this  process  before  a  firm  solution  was  ob- 
tained. 

The  Wiener-Hopf  method,  treated  by  Frankel  and  Goldberg  (3),  offers 
a  general  solution  to  the  one -dimensional  transport  problem  in  a  single 
serai-infinite  medium »   However,  it  is  available  only  as  an  exceedingly 
difficult  double  numerical  quadrature  which  could  not  be  handled  for  lack 
of  computing  facility.  A  solution  is  given  for  the  flux  at  the  black 
boundary  only,  requiring  only  one  numerical  quadrature.   By  equating  this 
solution  in  terms  of  the  parameters  to  the  experimental  value  of  the 
"contact"  flux,  unique  values  are  obtained  for  the  parameters.   These 
values  extrapolate  the  trend  of  the  nxjmerical  solutions,  and  are  in  that 
sense  to  be  regarded  as  supported  by  the  numerical  method. 

The  convergence  on  the  Wiener-Hopf  solution  may  be  illustrated  by 
the  successive  values  for  the  macroscopic  "total"  cross  section.  They 
ares   Numerical  -  2,9  cm~-'-j  2.5  cm"''"  and  2.3  cm""^;  Wiener-Hopf  -  2.10  crn'^. 


2„      SiJnmary 

Findings  may  be  tabulated  as  follows: 


Parameter  -  Water  at  24°  C. 

Values 

Diffusion  Length 

2.67  ±    .02  cm. 

Transport  Mean  Free  Path 

0.48  ±  .01  cm. 

Absorption  Mean  Path  (Equivalent) 

44.2  *  1.3  cm. 

Microscopic  Transport  Cross  Section 

62.4  =fc  1.5  barns/mol. 

Microscopic  Absorption  Cross 
Section  (Equivalent) 

0.68  ±  0.01  barns/mol. 

Ratio  -  Boundary  Flux  -r  Asymptotic 
(Exponential)  Flux 

0.186 

Ratio  ~  Boundary  Flux  -5-  Asymptotic 
(Hyperbolic)  Flux 

0.87 

Extrapolation  Distance  -  Xq 

0.32  cm. 

Ratio  -  Extrapolation  Distance/Total 
Mean  Path 

0.67 

Comparison  of  the  cross-section  with  the  energy-dependent  values  for 
protons  in  water  by  Bacher(l)  and  Melkonian  (7)  indicate  an  "average" 
neutron  energy  for  transport  of  perhaps  0.05  electron  volts  in  contrast 
with  the  Maxwell-Bolt zmann  mode  of  0.025  e.v.   This  is  not  surprising  since 
it  is  expected  that  transport  is  largely  effected  by  the  high  energy  neutrons - 
which  have  longer  free  paths.   In  addition  it  must  be  remembered  that  the 
neutron  flux  is  weighted  by  a  velocity  term  so  that  the  average  flux  corres- 
ponds to  a  higher  energy  than  does  the  average  neutron,  and  transport  weights 
the  cross-section  by  the  flux  rather  than  by  the  neutron  density.   The  l/v 
absorption  in  hydrogen  provides  a  preferential  sink  for  low  energy  neutrons, 
thus  the  thermal  distribution  is  displaced  from  the  Maxwell -Bolt zmann  toward 
higher  energy. 


CHAPTER  II 
TRANSPORT  THEORY  AND  METHODS  OF  ANALYSIS 

To  indicate  the  assumptions  involved  and  to  portray  the  physical  model 
which  was  analyzed,  it  may  be  well  to  review  the  fundamental  theory  of 
neutron  transport o   It  must  be  emphasized  that  the  model  of  neutron  trans- 
port for  neutrons  in  energy  equilibrium  with  the  medium,  at  energies  well 
below  the  atomic  binding,  differs  radically  from  that  for  energetic  neutrons 
(En  >  1  ev)  in  media  at  normal  temperature. 

In  the  latter  case  the  energy  transfer  in  neutron -nucleus  interactions 
is  uniformly  from  neutron  to  nucleus,  the  nuclei  interact  as  free  particles, 
and  since  the  atomic  vibrational  energy  is  comparatively  small,  the  nuclei 
may  be  regarded  as  initially  at  rest.   None  of  these  statements  hold  for 
the  former  case. 

In  the  present  problem.,  by  definition,  the  net  average  energy  transfer 
between  neutrons  and  medium  is  very  small  and  due  solely  to  the  preferential 
absorption  of  the  slower  neutrons.   The  protons  are  bound  in  the  molecule 
with  energy  larger  than  the  neutron  energy  and  the  nuclei  cannot  interact 
as  free  particles.   The  thermal  motion  of  the  atoms  and  of  the  systems  in 
which  they  are  bound  is  not  negligible,  and  so  affects  not  only  the  mechanics 
of  the  collisions  but  also  the  frequency  of  occurrence  since  the  collision 
probability  is  a  function  of  the  flux  of  nuclei  as  well  as  of  neutrons. 

Inelastic  interaction  of  neutrons  with  rotational  and  vibrational 
levels  of  the  water  molecule  are  probably  present.   The  rotational  levels 
are  quite  numerous  and  begin  at  wave  number  23.8  cm"-'-,  (about  .003  electron 
volts)  according  to  Wu  (9). 

Quantitive  discussion  of  the  scattering  mechanism  in  a  system  as 


complex  as  this  is  not  in  the  scope  of  this  paper.   However  to  avoid  mis- 
interpretation of  the  results,  it  is  necessary,  to  consider  at  least  the 
question  of  whether  the  scattering  is  symmetric,  that  is  whether  all  values 
of  the  angle  of  deflection  of  the  neutron  in  the  Laboratory  System  are 
equally  probable,  when   the  neutrons  are  in  energy  equilibrium  v/ith  the 
riedium. 

Referring  to  Fig.  (l),   i?   and  ^  are  the  velocities  of  neutron  and 
molecule  respectively  before  collision,  in  the  laboratory  system,  and  CJ 
the  angle  between  then:.   To  transfer  to  the  center  of  ffass  system  a  vector 
velocity  V  is  added  of  such  magnitude  and  direction  as  to  reduce  the  total 
momentum  to  zero.   Using  unity  as  neutron  mass  and  M  as  molecule  mass  we 
b^ve : 

Vo   +    V^    +    M(v  cos  cj    -I-  V^)    =0 


Vj    +   M(v   Sin  w    +  Vy)     =0 


V.  =  - 


M  I 


and  V  /  M 

Vy  =  -  ]^^prr  ^  ^"^^ 

where  \f      and  y   are  the  scalar  components  of  V  . 
X        ^ 

The  transformation  vector  is  thus  seen  to  resolve  into  two  vectors,  one  of 
magnitude  — — —  v     and  opposite  to   v    in  direction,  and  one  of  -^ag- 

M+l 

nitude        I  ^  and   op^josite   to      ^        in   direction,      The    vectors     ^' 

and   v'    are    .ho   velocities   before   collision   m  '^ ,-.  •    systei     and      v©        and  V 
are   -che   CM.   velocities   afler   collis-!on.      li"  the   scatter  ra      :^.    spherically 
symjTietric   ia  C.   m.   system,    all   values   of  cosy,    the   deflection  angle,    are 
e^"'^^  \y   i.  •^'^i^'l'' .      iC   c^  ^?jir    ',hc-  -neu-^^'o,^   veloGil.y  (^  ,     •^•z-t  collision,    in 

6 


Laboratory  system,,  we  add  (-V)  vectorially  to  %"  .     The  probability  of 
cj  is  sin  u)  ^  and  is  therefore  symmetric,  so  on  the  average  the  contribu- 


tion  of  ! V   inserts  the  direction  of  the  incoming  neutron  as  a 

component  of  u  .   Thus  the  scattered  neutrons  will  emerge  from  collisions 
preferentially  in  the  direction  of  their  pre-collision  motion  with  a  weight 
increasing  with  neutron  energy. 

The  preferentially  forward  scattering,  as  indicated,  correlates  rather 
poorly  the  direction  of  each  path  with  the  preceding  one,  less  well  with 
the  next  but  one  preceding,  and  after  a  large  number  of  intervening  scat- 
terings the  correlation  vanishes.   One  may  attempt  to  obtain  the  number  of 
neutrons  passing  through  a  small  volume  increment  (  dV)  by  summing  up  those 
neutrons  scattered  in  all  other  volume  increments  ( d  V )  in  the  space  which 
reach  the  first  increment  in  free  flight.  However,  since  the  number  of 
scattered  neutrons  leaving  ol  Y  in  the  direction  of  d  Y  depends  on  the  aver- 
age direction  of  the  neutrons  entering  dY  ,  this  average  direction  must  be 
found  by  another  summing  up  of  previously  scattered  neutrons  and  so  on. 
The  mathematical  analysis  of  such  a  situation  yields  infinite  order  multiple 
integrals,  so  simplification  is  necessary.  The  difficulty  does  not  occur  if 
the  scattering  is  isotropic,  so  an  obvious  suggestion  might  be  to  approxi- 
mate the  system  with  a  model  in  which  the  scattering  is  isotropic. 

If  one  considers  a  neutron  scattered  isotropically  at  a  point  in  a  non- 
absorbing  medium  and  asks  for  the  probability  of  finding  it  at  any  point  in 
the  space  after  the  lapse  of  a  long  time,  it  is  clear  that  the  resulting 
probability  field  is  symmetric  about  the  initial  point.  Analogously,  one 
considers  a  neutron  with  initial  motion  %     scattered  at  a  point  emisotrop- 
ically  in  the  Lab  system  but  isotropically  in  CM.  system. 


The  probability  of  finding  the  once  scattered  neutron  at  position  r 
is  then  symmetric  about  the  moving  center  of  mass  whose  velocity  is 

— — — —   V,     .   However,  the  probability  of  finding  the  second  scattered 
neutron  at  F   is  only  azimuthally  sj-mmetric  since  each  spherical  shell  of 
second  scattered  neutron  probability  is  immediately  again  retarded  by  the 
factor  ! ,  corresponding  to  the  motion  of  the  center  of  mass  of  the 

M  +  l 

second  scattering  systemic   Thus  the  second  scattered  neutrons  which  have 
made  long  flight  after  first  scattering  will  be  found  on  the  average  farther 
along  in  the  direction  of  V©  than  those  v;hich  make  short  flights.   In 
spite  of  the  forbidding  complexity  of  the  picture  some  information  can  be 
obtained  about  the  probability  distribution  after  a  large  number  of  collisions. 
The  path  lengths  X.     ,  the  azimuthal  scattering  angles  -li  ,  and  the  cosines 
of  the  polar  scattering  angles  u  ,  are  every  one  statistically  independent, 
and  may  be  averaged  independently.   The  centroid  of  the  distribution  is 
obviously  found  in  the  direction  of  V^    and  may  be  located  by  first  averaging 
over  the  azim:uthal  angles  starting  with  the  last  flight.   One  obtains 
(>f,/^,  +   ^2/^2/^1  -f  .  .  .  >?n  n, /^r>*^'")  •  Averaging  over  i^  and 
u^  one  has 

The  required  first  moment  is  then  )^  {  Tt<    +M^-j-/I^  +  ...  >u'^  +  ...) 

which  is  the  expansion  of    ^i/*-            ,   One  concludes  then  that  the  average 
position,  after  very  many  scatterings,  of  a  neutron  projected  from  a  point 
in  a  medium  is  at  a  distance    ^s('  "*"  t      ^  ^     °^    "i ^ —      from  the 

\  -M  I  -/A 

initial  point.   This  parameter  is  called  the  transport  mean  free  path.   The 
higher  moments  of  the  probability  function  are  not  available  in  general. 


F.  Joliot  (5)  has  corr.puted  the  probability  function  itself  for  the  n'^^ 

isotropic  scattering  in  a  quite  complicated  series  of  Herrnite  polynomials. 

One  may  now  reasonably  postulate  that  the  scattering  mechanism  in 

an  anistropically  scattering  medium  may  be  rather  well  reproduced  by  an 

isotropic  model,  provided  that  the  scattering  length  be  weighted  by  the 

factor   1    .   This  ignores  the  defect  that  the  higher  m.oments  are  not 

1  -JL 
properly  scalede 

It  must  be  emphasized  that  capture  has  not  been  considered  and  that 
the  model  is  not  correct  for  large  capture.   It  then  remains  to  be  demon- 
strated that  the  model  is  satisfactory  if  capture  is  small  but  not  negligible. 

The  simplified  transport  model  now  may  be  more  clearly  defined.   It 
consists  of  a  three  dimensional  array  of  fixed  lattice  points  at  which  the 
neutrons  scatter  isotropically,  or  are  absorbed. 

The  scattering  length  in  the  rrodel  is  the  transport  free  path  in  the 
real  medium  and  is  in  general  not  known.   The  absorption  free  path  in  the 
model  is  not  in  general  the  same  as  the  absorption  length  in  the  medium  but 
is  probably  nearly  the  same  if  absorption  is  small.   A  more  useful  pair 
of  parameters  are  the  inverses  of  the  free  paths  -  the  macroscopic  cross 
sections  -  defined  as  the  probability  per  centimeter  of  path  of  the  defined 
event.  To  prove  the  relation; 

~dl=  lid^  Where  I  is  the  beam  intensity 

Let        I^  =    I      ,         I   =    fc-^-^ 

The   mean  value   of    J(      is:  i  9  tT^^  Ai 

X  =  -^ 


z 


A  macroscopic  absorption  cross-sectior  is  similarly  defined. 
The  following  notation  is  to  be  used: 
n     -  Neutron  density,  cm''^ 

^    -  Neutron  scalar  flux  density,  called  simply  "flux,"  cm~^  sec"^ 
X  .   -  Mean  macroscopic  transport  cross  section,  cm" 
y  -   Mean  macroscopic  absorption  cross  section  cm" 

y    -  Mean  total  cross  section,  cm"'^ 
N    -  Lattice  point  density  of  the  medium,  cm~^ 
<r    -  J'^icroscopic  total  cross  section,  cm 
X  -   Distance,  cm. 

2     -  Distance  in  mean  paths 
r    -  Position  vector 

One  asks  for  the  number  of  neutrons  scattered  in  incremental  volume 
at  position  r  and  crossing  unit  normal  area  at  r'     per  second.   The  total 
scattered  at  F  is    ^(r)2]+olV   -  all  cross  sections  being  independent 

of  r   in  the  homogeneous  medium  -  and  of  these  d  w  leave  r  in  d  toi 

about  solid  angle  (a)  .   Since  the  area  subtended  at  r'     by  dcj  at  r   is 

(r'--   r)  du>    }  a  division  by  this  factor  is  indicated,  and  since  the 
beam  in  any  dw  is  attenuated  by  absorption  and  re-scattering  as 

^-Lx\r  -r\   ^-L.t^\r    -  r\         ^  multiplication  by  this  last  factor 

is  required.   The  required  number  is  then     ..        —     _2  r'  —  r  A'sr 

The  neutron  flux  at  r'    is  then  obtained  by  summing  up  the  neutrons 
scattered  at  all  r  and  crossing  unit  area  at  r       .   Thus 


In  the  problem  at  hand,  the  medium  is  considered  to  be  a  semi-infinite 
slab  -  infinite  in  horizontal  extent  and  downward,  its  upper  face  bounded 
by  a  plane  from  which  no  neutrons  return,  and  having  a  finite  neutron  current 
passing  upward  through  it.   The  coordinate  %     is  measured  from  the  "black" 
plane  and  normal  to  it.   In  this  system  the  flux  is  a  function  of  X  alone. 

One  requires  the  num.ber  of  neutrons  scattered  in  an  incremental  slab 
dX  at  X  and  passing  in  free  flight  through  unit  area  about  a  point  at  x'  , 
Referring  to  Fig.  2,  this  quantity  called  ^nx'-Xl)^"^   ^^   found  as  followsj 


Kda=     ^^    f  ^    (ZTrrdr 


dx) 


K      - 


rdr  =    pdp 


-'p=|)c-t| 


zK-x| 


The  integral  is  terrred   E,  (jj-x'-xj)       >  "the  exponential  integral, 
and  is  identical  with  the  function  — E{  (-x)   ^^  tabulated  by  Jahnke  &   Ehmde 
and  others. 

The  flux  at  x'  is  obtained  by  summing  the  contributions  of  all  increments 
dx  at  %    as  follows; 


¥^0c]    =    ^y  (^(xj  E,(2|^'«^,)dx 


space 

This  linear  integral  equation  is  perhaps  more  conveniently  handled  in  mean 
free  path  dimensions  as  follows: 


^(^')  =  ^/^ri^)  f.(„--„/^ 


For  the  configuration  space  described,  the  limits  of  integration  are  zero 
and  plus  infinity.   However  as   -£  —   —^   15  the  contribution  to  the 
integral  between  negative  infinity  ^nd  zero  becomes  negligibly  small,  qnd 
the  lim.its  may  be  taken,  with  vanishing  error,  over  all  values  of   E 
In  this  asymptotic  for-^^  the  equation  has  a  simple  hyperbolic  solution. 

The  usual  method  of  treatment,  to  indicate  the  form  of  the  asymptotic 
solution  (as  in  Frarkel  &  Nelson  (4)),  involves  transforming  the  origin  to 
write: 

and  expand   (^/g'^gj   in  Taylor's  series  about  H  J 

/oo 
-9^    C  _i  _        is  the  m"^^  m.oment  of  the  kernel 


about  its  origin,  and  this  vanishes  by  syn-^ietry  for  m   odd.  One  has  then 

neglecting  terir.s  above  the  third  order o  Rearranging,  one  has  the  one- 
dimensional  differential  equation  of  elen.entary  diffusion  theory  -  the 
Helmholtz  equations 


for  Twhich  the  solutions  are  hyperbolic  if    Til  q  ^  — = —   ,  circular  if 
Iff)   >   ^^      ,  and  linear  if   1110  =  -^ —        .   Now 

Tn  o  ~  /  E|/   vdZ  =  2    5  and  if  absorption  is  present,  then  Z^^t  » 

so   ~Y\  -  9  ^     -  '      Thus  we  expect  hyperbolic  solutions,  so  a 

(Mo  -  ^  ^  y 

^     —  ^^ 

trial  function  U)f  \  =  ^  is  substituted  in  the  integral  equation,  and 

the  integration  perforir:ed  at  the  point  Z    —  O 

z,    r"^  jf2 


It 

2Z 
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~      T        I   cosh  (^?j  E, /g)  dZ  which   is    integrated 

by  parts   after   putting     E|  (^)    in  integral    foriTo 


I 


COSh(c^g;d&    /        ^  d 


The  first  integral  vanishes  f or  j?<  |  since  at  upper  lin;it  ^1(9)      is  of 

order  k,  and  at  lower  lirdt  — '-^ —      goes  as  *  and  the  singularity  of 

E  I  (•5^   is  logarithmic ,   Writing  the  second  integral  in  exponential  form 

i,    z^   A°°,-;(i-j()-,-zO-^x)  d2 
2*1     /  z 


and  separating  the  range  of  integration  into  intervals: 


■^fl^.  -f^...£4^.-f»j 


the  third  and  fourth  integrals  are  seen  to  be  each  Ej^Cl)  and  cancel,  thus 

I  -    It     o    LtL 


2:   ■     ^Xv,i±^      t..K-'je 


expanding  in  power  series 


^=      I    ^41+4:+...+  Jl!l 


5  5       '   •••  2h+l 


Inverting   the    series    in    v(       we   obtain 


a2 


■laf  .o.,  /'I*'\* 


function   in   crrio    meas 


an  eigen  function   in  irenn   free 
ure   is      fe  where    k  =  «>  ^     . 


Note   that     e"*  is   an  eigen  function   in  irenn   free   path  measure   and   the 


Given  an  experirrer.tal  value  of  K  (which  is  siirply  the  logarithmic 
slope  of  the  flux  distribution  plotted  against  "X  )  one  may  choose  a  value 
of  2.  ,  compute  ^  byjv=/^,  and  using  the  series  in  J^   (or  any  conveniently 
obtained  solution)  compute  ^y^     °   Returning  to  the  integral  equation 


^f^)  =  #-/W;E,(,,,_^,^di' 


one  may  insert  the  computed  value  of  /^c-  •   It  is  now  required  to  com- 

pute the  eigenfunction  ^H\   corresponding  to  the  selected  value  of  ^  , 
A  choice  of  method  is  available. 

LeCaine  (6)  has  treated  the  problem  by  a  variational  method,  presenting 
tabulated  values  for  the  flux  normalized  to  unit  flux  at  the  black  plane 
for  various  values  of  ^^y  «  A  rather  rough  solution  for  the  parameters 

is  available  by  using  these  results,  and  indicates  that  2m     lies  between 
2oO  cm~^  and  2o4  cm'-^. 

It  was  then  decided  to  attempt  solution  by  a  direct  numerical  method, 
which  was  carried  out  to  completion.   In  theory  the  method  goes  as  follows: 
Given  the  integral  equation  as  before 

%,  =^/^(2')  £,(,,,.„)  di' 

The  integral  is  replaced  by  a  Riemann  sum  equation 
which  may  be  put  in  matrix  notation 


15 


It 

11 


■A2 


too  tg,  t,j2  ~  "■ 

E,c  E„   E„ 

^20     ^21       "^22       -        -        - 


-   El 


where     bj_^    is    the   absolute  value   of  the   exponential   integral   of  argument 
(i.—  n)  A  E         •      A't   ■the   argument   i   =  n  the   exponential   integral    is    singular, 

so  a  raean  value  was   computed  for   this    range.      One   requires: 


4f 
"   2 


-^      j-<».^i 


h 


The  integrating  scheme  is  the  equivalent  of  the  trapezoidal  rule  but 
employing  the  function  of  r^ean  argu-nent  instead  of  the  chordal  mean  as 
the  suirmaticn  eleirent, 

On.e   proceeds  by  inserting  a  tabular  trial  fimction,  perforring  the 
matrix  operation  as  indicated  and  emerging  -./ith  a  new  and  bettered  function 
which  is  iterated  until  the  variation  prodiiced  by  iteration  is  insignificant. 
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In  practice  ninor  departures  from  the  straij^ht  forward  process  were  intro- 
duced to  reduce  the  labor.   The  principal  one  of  these  consisted  in  iterating 
for  only  about  10  terms  of  the  flux  vector  and  writing  the  cojTiplete  vector 
by  smoothing  the  first  differences. 

A  variant  of  this  scheme  was  considered,  which  involved  writing  the 
equation  in  a  perturbation  forms 

where  ifo^   is  the  asymptotic  solution,  8  is  the  boundary  depression,  and  c 
is  the  perturbation  of  the  operator  Ej  required  to  extend  the  integration 


over  all     2      . 


hi)  =  ^iirj  { [E'(i?'-zi)  ~ ^(i^'-zi)] «(?)  +  e V.X I dt 

Tl 


or   in  Rieraann  sum  form 


K        21 

The   equivalent  matrices,    all    ini'inite   order  are; 


k- 


i- 


-l-Z    ^-2-1    ^.l\  ^-11 

E-i-1  E.-i-»   t.|,  t.|2  - 

E...  E,.,    0  0  0 

^2-2   K^    0  0  0 

-      -     0  0     0 
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E-£ 


0 

0 

0 

0 

0     0     0 

0 

0 

0 

0 

0 

0      0     0 

0 

0 

0 

0 

0 

0     0     0 

0 

0 

0 

0 

0 

0     0     0 

0 

0 
0 
0 
0 

0 
0 
0 

0 

0 
0 
0 
0 

- 

0 

0 

0 

- 

-      -     _ 

- 

The  first  method  has  the  disadvantage  that  the  term  products  converge 
rather  slowly  for  increasing  2  so  that  each  term  of  the  iterated  vector 
represents  very  irany  term  products  and  the  computation  is  quite  laborious. 
The  second  method  avoids  this  trouble  since  S   vanishes  strongly  with  in- 
creasing 2  ,  but  S    varies  so  rapidly  near  (i=  O)  that  trapezoidal  inte- 
gration is  valid  only  for  a  very  small  interval  Ab  .   On  balance,  method 
one  appeared  most  feasible  and  was  choseu.  A  value  of  2.9  cm   was  obtained 
for  2.  with  a  very  fine  fit  of  the  solution  to  the  experimental  data. 

To  investigate  the  validity  of  the  integration,  the  process  was  then 
repeated  with  an  improved  kernel.   In  this  case  the  matrix  elements  were 
obtained  by  computing  mean  values  for  the  exponential  integral  with  a  unit 
weight  function  for  each  Az^o   This  process  amounts  to  correcting  for 
variation  of  the  kernel  over  the  interval  but  ignoring  the  variation  of  the 
function.   A  value  of  Z  =  2.5  cm"-^  resulted. 

A  further  improvement  in  the  kernel  was  then  tried.   One  assumes  that 
the  function  nay  be  linearized  with  small  error  over  the  interval  Ai 
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The  matrix  equation  is  then  written; 
rA2 


where  2^  is  the  lower  boundary  of  the  n'^^  interval,  Az  is  the  width  of 

the  interval,  A^n  is  the  difference       ^(j.       \     —  C^/z  \\  »    ^^^ 

^     is  the  difference   \i   ~  2^    •   The  equation  nay  be  separated  into 


two  terms  as  follows; 


i\)'' 


The  first  integral  (letting  y  =  |2c-2„  +  ^(  =  |Kr,H-f|) 

The  second  integral,  dropping  absolute  values: 


(t|-K„)E„„^ol 


=  4-[e-'^"(ltK„)  -e-''"'(|tK„J+CE,(,„^,)-K:E,(^j-K,I,,„ 

I-     constitutes   the    terms    of   an   operator  which  operates    on  the  function   iyn 
and  l^  an   operator   vrfiich  operates   on  the    relative    increment   of   the   function, 
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The  process  yielded  again  an  excellent  fit  to  the  data  but  for  a 
value  Z  s  2.3  cm"'^.   Since  the  computing  had  become  impractically  tedious, 
further  v;ork  on  this  line  v/as  dropped. 

Frankel  and  Goldberg  (3)  list  a  solution  for  the  boundary  flux  as 
follows  5 


%,  -  VHiuSr^) 


where  the  asymptotic  flux  goes  as  the  hyperbolic  sine  vanishing  at  an 
"extrapolated"  boundary  thuss 

The  extrapolation  distance  is  given  by: 

The  characteristic  equation  of  the  systeir  relates  i^   to    t   ,  so  th 
boundary  flux  may  be  expressed  as  a  function  of  it  alone; 


L{\-e) 


3        "^        15        "^      35       ■^"•(2n+)X2n-0         '      TTJif^ 
'^       ^        5       ^       7        ^  •  •  •   2n+l  I  +  2  (/?(?, 
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The  procedure  is  as  follows: 

1.  The  extrapolation  distance  is  con-puted  by  Simpson's  Rule  for 
a  few  values  of  ^  in  the  neighborhood  of  0.15  -  0.20  as  in- 
dicated by  the  numerical  solution,  and  plotted  against  X  , 

2.  Selecting  a  value  of  ^  and  the  corresponding  Z©   ,  the 
experimental  boundary  flux  is  normalized  to  go  as 

sinh  i?  (^  +  Bo)  at  large  argument.   This  new  normali- 
zation is  obtained  by  halving  that  flux  which  was  normalized 
to  go  asymptotically  as  fc     and  multiplying  by  fc   °   . 
The  result  is  (fp)   . 

3.  Insert  the  Y(p)    in  the  series  equation  and  solve  by  trial 
for  ,)^  , 

4.  Using  the  con-puted  values  of  iv  ,  return  to  stage  (2)  and 
iterate  the  process.   Convergence  is  rapid. 

The  described  process  yielded  L  =  2.10  ±.  0.05  cm"^ . 


CHAPTER  III 
EXPERHIENTAL  PROCEDURE 

The  medium  in  vi^hich  the  experiment  was  conducted  was  the  approxi- 
mately 6'  X  6'  X  6'  water  tank  which  is  ir.ounted  on  the  thermal  column  of 
the  Oak  Ridge  Graphite  Reactor,   A  current  of  thermal  neutrons  passes 
upward  through  the  tank  after  filtering  through  the  graphite  column  below. 
The  water  content  is  ion-exchange  demineralized,  and  its  nuclear  properties 
are  substantially  those  of  distilled  water. 

A  disk  of  aluminum,  faced  with  ,020  inch  cadmium,  and  approximately 
10"  in  diameter,  carried  on  aluminum  legs  about  10"  long,  was  lowered  with 
a  rod  to  a  reproducible  central  location  on  the  tank  bottom.  A  foil  holder 
of  aluminum  and  lucite,  fastened  to  the  cadmixom,  face  was  provided  with 
adjustable  mountings  for  positioning  foils  along  the  central  normal  to  the 
plane.  A  foil  was  also  cemented  to  the  center  of  the  cadmium  plane  for 
each  exposure.   The  foil  mountings  were  of  lucite  of  minimum  section  con- 
sistent with  stiffness,  and  a  film,  of  paraffin  v/as  employed  as  an  adhesive 
to  mount  the  foils,  thus  introducing  no  non-hydrogenous  material  in  the 
immediate  vicinity  of  the  foils.   The  equipment  is  shoivn  as  Fig,  (3), 

Indium  foils  were  used,  weighing  25  milli-grams  nominal,  cut  from 
50  mg/cmi  stock.   These  were  cleaned  of  the  adhesive  wax  after  each  use, 
weighed  and  allowed  to  decay  for  at  least  48  hours  before  re-use. 

Foil  counting  was  done  with  an  end  mica-window  Geiger-Mueller  counter, 
connected  through  a  plate-loaded  pre-amplif ier  stage  (of  gain  -^   12)  to  a 
scale  of  64.   Counter  and  pre-amplif ier  were  enclosed  in  a  2"  lead  shield, 
which  reduced  background  to  a  negligible  fraction  of  the  foil  activity. 
Sufficient  counts  were  obtained  in  each  counting  run  to  reduce  statistical 
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fluctuation  to  a  negligible  level,   A  nunber  of  determinations  of  counter 
resolving  time  were  made  by  following  the  decay  of  54  min  indium  116  over 
about  5  hours o   An  average  value  of  240  yuseca  was  obtained  by  solving 
the  following; 

A.      ,      Aa 

where  A  :;  extrapolated  activity 

N  =  mean  counting  rate,  sec~^ 
0"=  resolving  time,  sec. 

A  correction  factor  l/l-N3"   was  plotted  against  counting  rate  and 
applied  to  all  counting  runs  to  correct  for  coincidence  losses. 

For  exposure,  the  foil  holders  (after  mounting  the  foils)  were  adjusted 
approximately  to  position  with  respect  to  the  cadmium  plane  using  a  common 
scale.  A  vernier  height  gauge  was  then  used  to  determine  the  exact  foil 
position  (estimated  probable  error  ±.    ,05  mm.).   The  plane  and  foil  holder 
assembly  was  then  lowered  quickly  to  the  tank  bottom  and  exposed  for  times 
ranging  from  5  to  15  minutes  depending  on  foil  position  to  obtain  reasonable 
counting  rates.   No  attempt  was  made  to  measure  either  exposure  or  foil 
"cooling"  times  accurately  since  one  foil  always  and  two  foils  usually  were 
exposed  at  fixed  positions  to  monitor  the  activity  of  the  variable  foil. 
The  positions  at  x  -  0  and  xS  5.18  cm.  are  the  monitor  positions  and  data 
at  these  points  represents  an  average  of  about  15  measurements.   In  addition 
to  the  exposures  made  with  axial  foils,  two  exposures  were  made  with  foils 
set  in  a  plane  parallel  to  the  cadmium  to  demonstrate  that  the  geometry  was 
actually  one  dimensional.  A  plot  of  the  diair.etral  variation  of  flux  appears 
as  Figo  (6).   The  maximuir.  deviation  from  the  axial  flux  amounts  to  about 
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5%   reduction  at  8  cm,  radius.   This  is  considered  satisfactory.   To  test 
the  flux  for  fast  components,  a  cadmium  difference  on  Indium  was  taken. 
The  ratio  of  activity  in  Cd,  to  activity  in  the  "shadow"  of  the  cadmium 
envelope  (a  most  conservative  condition)  exceeded  5,000  -  therefore  foil 
activation  due  to  the  1,44  volt  resonance  may  be  safely  neglected. 

To  analyze  the  data  by  the  numerical  method,  it  vra.s  first  plotted 
against  the  coordinate  x  in  cms.  on  semi-log  paper  and  the  fair  curve 
drawn  by  eyeo  A  clearly  exponential  branch  was  obtained  beginning  at 
X  «-^4.8  cmo   The  decay  constant  k  of  the  exponential  was  solved  by 

a.  —  b 

and  fixed  at  k  s  0.375.   The  scale  of  the  values  of  (^(x)  was  then  adjusted 
so  that  (P   (asyir^ptotic)  =   ^  *    .A  trial  value  of  X   was  then  selected, 
and  the  data  plotted,  faired  and  tabulated  as  a  function  of  2  in  mean  path 
measure.   The  eigenvalue    ^  -.   was  computed,  the  ^if\     A-ai^    tabulated 
on  a  strip  for  Az  =  "S,  and  numerical  iteration  of  the  equation  carried 
out  until  it  was  obvious  that  the  solution  near  the  plane  at  E  s  0  would 
not  fit  the  data.   The  trial  was  then  terminated  and  another  trial  made  for 
a  different  value  of  Z.  , 

As  the  solutions  were  improved  it  becajne  clear  that  a  correction  was 
necessary  for  the  remainder  contributed  by  the  ter.ns  beyond  the  last  tabular 
value  of  the  Ei .   The  correction  was  inserted  by  computing  the  remainder 
with  respect  to  a  proper  exponential  function  and  dividing  by  the  value  of 
this  function  at  the  first  untabulated  term  of  the  E^^ .   This  number  is 
then  used  to  tenninate  the  E-^   tabulation.   Since  the  last  tabulated  value 
of  the  El   was  for  inAZ  =7.4  one  obtains; 
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■'-ic: 


"^^m'^^ 


which  may  be  integrated  by  writing  the  Y.^   as  an  integral,  then  inverting 
the  order  of  integration  to  yield 


which  must  be  computed  individually  for  each  value  of  u.    , 

Iteration  was  carried  out  by  computing  only  every  fifth  term  of  the 
trial  function J  then  interpolating  the  new  function  by  making  the  first 
differences  regular o   The  final  integration  was  done  in  detail  in  the 
region  close  to  the  black  boundary. 

Plots  of  the  flux  as  a  function  of  distance  from  the  black  boundary 
appear  as  Figo   4)  and  (5)o   The  curves  were  obtained  by  using  the  Le  Caine 
results  to  interpolate  the  data  near  the  boundary^  and  the  numerically 
iterated  function  in  the  region  more  than  0«5  cm,  from  the  boundary.   The 
latter  process  is  considered  reasonable  since  the  solutions  in  this  region 
were  invariant  to  changes  in  the  integrating  scheme,  they  are  constrained 
by  the  approach  to  the  asymptotic  function,  and  they  fit  the  data  quite  well. 
Desirable  extensions  to  the  analysis  might  include: 

(1)  Computation  of  the  'Tiener-Hopf  solution  for  perhaps  15  or 
20  points  to  firm  up  confidence  in  the  fit. 

(2)  Extensive  analysis  of  the  numerical  iteration  method  to  estab- 
lish the  level  ^f  computing  cte "tail  necessary  to  yield  accurate 
solutions  using  the  true  values  of  the  parameters. 

Both  propositions  require  computing  machine  service  to  reduce  program  time 
to  a  practical  value, 

25 


i«inii'S'" 


R6.I 


VELOCITY  DIAGRAM  FOR  TWO  BODY  COLLISION 
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FIG.    3 
FOIL  MOUNT  AND  CADMIUM  PLANE  ASSEMBLY 
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APPENDIX    II 


The   flux   is    tabulated  as   a   funcbion  of    ^  for  tv\ro  values    of   the 
total   cross    section.      These   are   solutions   corresponding    (for  <^  r   2,3) 
to  the  most   elaborate  numerical   treatment,    and    (for  ^r   2.1)   to  the  best 

reconciliation  of  i'7iener-Hopf ,    variational,    and  numerical    and  experi- 
mental  solutions. 


^ 

^(1=   2.1) 

(j^il:.   2.3) 

2 

^(5=   2.1) 

^(f  =   2.3) 

0 

1.86 

1.86 

4.2 

18.00 

16.45 

0.2 

2.71 

2.70 

4.4 

18.90 

17.20 

0.4 

3.44 

3.45 

4.6 

19.80 

18.00 

0.6 

4.15 

4.14 

4.8 

20.75 

18.80 

0.8 

4.84 

4.80 

5.0 

21.75 

19.70 

1.0 

5.52 

5.44 

5.2 

22.75 

20.60 

1.2 

6.20 

6.07 

5.4 

23.80 

21.50 

1.4 

6.88 

6.69 

5.6 

24.90 

22,40 

1.6 

7.57 

7.30 

5.8 

26.00 

23.25 

1.8 

8.29 

7.95 

6.0 

27.10 

24.20 

2.0 

9.04 

8.59 

6.2 

28.20 

25.20 

2.2 

9.84 

9.25 

6.4 

29.50 

26.20 

2.4 

10.64 

9.93 

6.6 

30.70 

27.25 

2.6 

11.49 

10.60 

6.8 

32.00 

28.30 

2.8 

12.30 

11.30 

7.0 

33.40 

29.45 
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7.2 

34.80 

30.55 

3,2 

13.85 
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7.4 

36.10 

31.75 

3.4 

14.65 

13,50 

7.6 

37.60 
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15o45 

14.25 

7.3 

39.10 

34.20 

3,3 

16.30 

14.95 

8.0 

40.80 
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17.15 

15.70 

38 


SI 


r 


15537 
1 Thesis 
S5*    ^'""Transport  parameters 
'        for  thermal  neutrons 

in  water. 


